PIECEWISE CERTIFICATES OF POSITIVITY FOR MATRIX 

POLYNOMIALS 
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Abstract. We show that any symmetric positive definite homogeneous ma- 
trix polynomial M G R[a;i, . . . , Xn]™'^"' admits a piecewise semi-certificate, 
i.e. a collection of identites M{x) = fij(x)Uij(x)'^Uij{x) where Uij{x) 
is a matrix polynomial and fij{x) is a non negative polynomial on a semi- 
algebraic subset Si, where R" = U^_-^5'i. This result generalizes to the setting 
of biforms. 

Some examples of certificates are given and among others, we study a vari- 
ation around the Choi counterexample of a positive semi-definite biquadratic 
form which is not a sum of squares. As a byproduct we give a representation 
of the famous non negative sum of squares polynomial x'^z-^ + z'^y^ + y'^x^ — 
Zx^y'^z^ as the determinant of a positive semi-definite quadratic matrix poly- 
nomial. 



1. Introduction 

The Hilbert 17-th problem asks if a non negative polynomial fix) € . . . , a:„] 

is a sum of squares. The answer given by Artin provides an identity with denomi- 
nators (or rational functions), namely there are two sums of squares of polynomials 
q(x) and p(a:) such that q{x)f{x) =p{x). 

Such an identity is a called Positivestellensatz and can be seen as a certificate 
which algebraically proves the positivity of the polynomial f{x). For a general 
polynomial, the denominator q(x) is necessary. 

One may also consider a "relative" version, namely when f{x) is non negative on 
a basic semi-algebraic subset S — {x £ M.^ \ fi{x) > 0, . . . , fr(x) > 0}. This search 
of certificates received a lot of contributions especially for those semi-algebraic 
subsets S where the Positivestellensatz is denominator free. Essentially it happens 
to be possible when f{x) is positive and when the description of S satisfies some 
archimedean property (which can be seen as a "strong" compactedness assumption). 
Among all the available references, we mention the seminal works of Schmiidgen 
P] and Putinar 0. 

Another natural and often difficult question is, when it is possible to obtain a 
certificate of positivity, to study its complexity, i.e. the number of squares that are 
needed. 

In this paper, we are mainly interested in matrix polynomials, namely polyno- 
mials whose coefficients are matrices. Let x — (xi, . . . , a;„) and R[xY 
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of all square matrix polynomials of size m with entries in M.[x]. In that context, 
by a square we mean a expression of the form U{x)^U{x) where U{x) e R[a;]™^™, 
sometimes we will also use the terminology of hermitian square. 

Given a symmetric matrix polynomial M{x), the natural question is the following 

Assume that M{x) is positive semi-definite {psd in short) for all x e M". Is M{x) 
a sum of hermitian squares ? Or in otherwords, do we have a positivity certificate 
for M{x) ? 

It is well known that, when n — 1, then M{x) is a sum of two squares (see for 
instance [4 for a proof), i.e. M{x) = U (x)'^ U {x) + V {x)"^ V (x) where U{x),V{x) S 

There is also another case where we have a positive answer : the case of bi- 
quadratic forms in dimension 2. Namely, if m = 2 and if the entries are homoge- 
neous polynomials of degree 2, then M{x) is also a sum of squares. 

In fact this last case can be seen as a particular case of the homogenized version 
of the previous one (see [2] for this fact and a proof). 

Unfortunaltely these are essentially the only cases where such a certificate exists 
for a matrix polynomial. For instance, it becames false even for biquadratic forms 
in 3 variables and dimension 3 (confer the Choi counterexample in [T]). 

We may also look at the relative situation : if M{x) is positive definite for all 
a; in a semi-algebaric subset S C M" whose description satisfies an archimedean 
property, then we get a similar certificate as in the polynomial case. This has been 
proved by Hoi and Scherer in [6] ; see also [7] for a different and more algebraic 
proof. 

In this paper, we are searching for certificates of positivity without denominators 
for matrix polynomials. We first study certificates which are sums of squares of 
matrix polynomial weighted by a psd (scalar) polynomial. More precisely, given a 
psd matrix polynomial M{x), we wonder if M{x) can be written as a finite sum 

(1) M(x) =^/,(a:)(7,(a;)^C/,(x) 

i 

where Ui{x) G M[x]™^'" and Ji{x) is non negative polynomial on M". 

We call it a semi- certificate since, roughly speaking, we have a certificate which 
is a hermitian sum of squares weighted by some coefficients which are non negative 
polynomials. 

Of course such a certificate may gives some others which will be totally alge- 
braic by using usual certificates for the psd polynomials fi{x). Beware that using 
denominators at this point could break any interest since we may readily get one 
such certificate with denominators simply by using the Gauss algorithm (indeed, 
it suffices to inverte the m — 1 first principal minors of the given matrix poly- 
nomial). Although, whene the matrix polynomial is positive definite, a relative 
version of semi-certificate with respect to a semi-algebraic subset S whose descip- 
tion is archimedean would produce a usual denominator free algebraic certificate for 
the matrix polynomial on S, using for instance Schmiidgen or Putinar certificates 
for all the fijS- 

Unfortunately such a semi-certificate as in ([1]) does not exist in general, even for 
a positive definite matrix polynomial, as it is illustrated in section 3. 
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Nevertheless, noticing that one may obtain such certificate locahy, we introduce 
the notion of piecewise semi-certificate, i.e. a collection of semi-certificates as in ([!]) 
with respect to a semi- algebraic covering of R". 

We show in section 4, that such a certificate exists for any positive definite 
homogeneous matrix polynomial. We show also that we may translate the notion 
of semi-certificate and our related results to the case of biforms. 

Some examples and counterexamples are given and for instance, in section 5, we 
study some variations around the Choi counterexample given in [T]. 

As a byproduct we give a representation of the famous psd non sum of squares 
polynomial x^z'^ + z^y^ -\- y'^x^ — 3x^y^z^ as the determinant of a psd quadratic 
matrix polynomial. 

Of course, a lot a work have to be done to better understand the framework 
of semi-certificate for psd non definite matrix polynomials. Moreover, concerning 
complexity, one may be interested in estimating the lenght of the hermitian sum 
of squares representation, and also in the number of pieces of the semi-certificate. 
Since we use a compactedness argument, no bound on the number of pieces of the 
semi- algebraic covering can be derived by our method. 

2. Notations 

2.1. Matrix polynomials and biforms. In all the paper, we consider a set of 
variables x = (xi, . . . , a;„) and the polynomial ring M\x\ whose associated field of 
rational functions is R(2;). 

A matrix polynomial A{x) is just a matrix whose entries are polynomials. For 
instance, for a square matrix polynomial of size m in the variables x, we write 
A{x) e K[cc]'"^™. 

A form of type (n, d) will denote an homogeneous polynomial of degree d in n 
variables. A biform of type (ni, di; ri2, ^2) will denote a homogeneous polynomial 
in ni + n2 variables which is of degree di with respect to a set of ni variables and 
of degree ^2 with respect to the other set of d2 variables. 

To a matrix polynomial A{x) = {ai,j{x)) e R[x]™^™ whose entries are all forms 
of type (n, d), we may canonically associate a biform of type (ni, di; to, 2) : 

fA{x,y)= ^ a,j{x)y^yj. 

l<i,j<m 

We say that a biform of type (ni, di; 712, ^2) is positive definite if f{x, y) > for 
all {x,y) G 5"^"^ X S*"^"^, where S^~^ denotes the ki dimentional unit sphere in 
R'^. Equivalentely, we may alternatively consider the product of projective spaces 
instead of unit spheres. 

If A{x) is a symmetric matrix polynomial in R[a;]™^'" whose entries are all forms 
of type (n, d), we say that A{x) is positive definite if the associated biform /^^(a;, y) 
is positive definite. We may also note that A(x) is psd if and only if fA{x, y) is psd 
on the product of unit spheres. 

2.2. Sum of squares and Hilbert 17th problem. Note that sums of hermitian 
squares when dealing with matrix polynomials corresponds to usual sum of squares 
in the terminology of biforms. 
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Let A{x) be a symmetric matrix polynomial in R[a;]'"^™ which is positive semi- 
definite for all substitution of x e M" . Then, the Gauss reduction algorithm gives a 
solution with denominators to the Hilbert 17-th Problem for a generic matrix poly- 
nomial : there exist r e N, Ui{x), . . . , Ur{x) e and fi{x) . . . , fr{x), f{x) G 
M.[x] psd polynomials such that 

r 

(2) f{x)A{x) - h{^)M^fUdx) 

i=i 

where f{x) can be choosen as the product of the m — 1 first minors of A{x) when 
they do not identically vanish. Moreover, by multiplying the identity ([2]) by the 
denomiators appearing in a positivity certificate for the /i's, we may assume that 
the polynomials fi and / are sum of squares. 

In our paper, we are mainly interested in positivity certificates without denom- 
inators, namely when we may choose / = 1. Hence, we first look for matrix 
polynomials which can be written as in ([T]) : 

r 

A(a:) =^/,(a;);7,(a;f ;7,(x) 

i=l 

where the fiS are non- negative polynomials (we cannot ask for them to be sum of 
squares since denominators are already needed for polynomials). 

One may want to measure the complexity of the certificate given in ([T]), simply 
by counting the number of squares : r. But, one may also prefer another way of 
counting. 

Noticing that C/^C/ = U'^{Ei+. . .+Em)U = {EiUY{EiU)+. . .+{EmU)^ (EmU) 
where Ei is the diagonal matrix whose only non null entry is 1 at the i-th posi- 
tion onto the diagonial, we may prefer to count "rank-one" sums of squares, i.e. 
when U{x) has rank one in ]R(a;)'"^'". The advantage of that count is that it coin- 
cides with the number of needed squares when we view our matrix polynomial as 
a biform. 

3. Semi-certificates 
The following examples shows that it is not possible in general to hope for ([T]). 

3.1. Some examples. 
Example 3.1. Let : 

1 + xy \ 

xy x"^ +y'^ J 

Of course M is psd for all (x, y) e M^. 

To show that M cannot be written as in ((1}, we may proceed by applying to 
the biform /m which is canonically associated to M) the technology of cages and 
Gram matrices developped in . But we prefer to produce here an elementary and 
self-contained argument. 

Assume that M = J2l=i hUfU^, where G {m.[x,y]f'^ and > on aU R^. 
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By decomposing into rank one sum of squares, we may assume that 

UfUi = ^ J^'^ °^2* ^ I where ai,6i e 

Identifying the (1, l)-entry, we get X]i=i /j^f — ^ + ■ Thus, fi does not depend 
on y : fi & R[x]. Since fi > for all x e R, we deduce that fi is a sum of two 
squares in R[x]. Thus, if we increase the integer r and if we change the UiS, we 
may assume that fi — 1. Identifying the entries, we get 

= 1 + ^^ 

By a simple degree consideration, the polynomial bi has necessarily the form 

bi = aiX + fi^y"^, with ai,^^ £ R. 

Likewise, the polynomial has necessarily the form 

ai = Ti + S^x with 7,^, Si e R. 

It follows a contradiction with the identity '^^^i ~ ^V- 

The previous pathology is not due to the fact that the matrix M is only psd, as 
shown by considering the following variation : 

Example 3.2. Let 

1 + x'^ + e{x* + y*) xy 

xy e(l + x") + + 

where e > 0. This is a positive definite matrix polynomial on all R^ (and even 
at infinity, i.e. the associated homogeneized matrix polynomial remains positive 
definite). 

Assume that N can be written as in ([T]) : 



iV = y^f^ix,y) 



Ai AB,, 



1^ A,B, B? 

r 



RkSk Si 



where the fi's are polynomials of degree at most 4, the gj,Rk, S^s are polyno- 
mials of degree at most 2, the Uj,Vj's are polynomials of degree at most 1 and 
the Ai,Bi^s are constant. Moreover, the polynomials fi,gj are assumed to be non 
negative on R^. Since their degrees are at most 4 and the number of variables is 2, 
we know that they are sum of squares of polynomials. Hence, we may assume that 

(3) N = y^( 



where 

Ri = ftj + biX + CjX^ + diy + Ciy'^ + fixy, 
5*1 = Q!i + PiX + ^ix"^ + 5.iy + Viy'^ + Hixy, 
{ai,bi,Ci,di,ei,fi,a^, Pi, -fi, 6^,1^1, Hi) G R^^ 
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Let us consider the following vectors in the usual euclidien space : 

a = (ai),5_= {h),c= {ci),d^ {di),e^ {ei),f ^ (fi), 
a = {a,),P = (A), 7 = = {5i),v = {i^i),fl = {^li). 

Identifying the non diagonal entries in ([3]), we get 

(4) h-S + P-d + d-fi + J-a^l 

where • denotes the usual inner product on R''. 

By identifying the diagonal entries in ^ and by the Cauchy-Schwarz inequality, 
we get 





= 1 - 2a • c 


< 


1 +2%/s2\/g2 


< 


p 


= -lOL-V 


< 




< 




= 1 - 2a • 7 


< 


1 + 2e 




d^ 


= —2 a • c 


< 


2V^ 






= 1 










= -27-P 


< 


2V^ 




P 


= -2c-e 


< 


2%/^\/i2 


< 


0? 


= e 









2Jl 



< 2e 



Again, by the Cauchy-Schwarz inequality, we have 



b-S < v /2(l + 2Vi)v ^ 
/3-J < v/2(TT2^ 
a ■ fi < \/Ve 
f-a < V2e 

Thus, if we take e small enough, then we get a contradiction to Namely, N 
does not admit a certificate as in ([1]). 

3.2. Semi-algebraic covering. Thus, even for positive definite matrix polynomi- 
als, a certificate ^ does not exist in general. Although, one may remark that such 
a certificate exists locally in our examples. 



Look at example 13. II When \x\ > \y\, then consider the identity 



xy 

„2 I „,4 



X xy 
xy y2 



1 + x 
xy X' + y 

whereas when \y\ > \x\, then consider 

1 + xy \ ( \ 

xy 



1 




+ (x" -y^ + y4) 



y 



xy 

9 9 

xy X y 




1 




1 



Hence we are going to change our definition by considering a semi-algebraic 
covering of our space. 

From now on, we will focus on forms, namely homogeneous polynomials. 

Definition 3.3. Let f{x,y) be a biform of type (n, c?i; m, c?2)- A semi- certificate of 
positivity is the data of a semi-algebraic covering R" — SiU . . .U Sr such that 



^x € Si, f{x,y) 



PIECEWISE CERTIFICATES OF POSITIVITY FOR MATRIX POLYNOMIALS 



7 



where each ,fi,j{x) is a form of degree ei.j which is non negative on Si and gf j{x, y) 
is a hiform of type (n, di — Ci^j] 771,^2). 

Roughly speaking, the semi-certificate is a piecewise identity which is a sum of 
squares with respect to one set of variables y and only psd with respect to the other 
set of variables x. 

For convenience, we may formulate what happens when we are dealing with 
matrix polynomials. The symmetric psd matrix polynomial M{x) S R[x]™^™ 
admits a semi-certificate if there exists a (finite) semi-algebraic covering (S'i)i<i<r 
of M" and 

(Vz)(3r, e nm e {1, . . . ,sj)(3/,,,(a;) e R[x]){3U^,,{x) e M[x]"X") 
such that 

(5) ^Af (x) = J2 f^A^)U'[Ji^)U^A^)^ and {Vx e 5„ f,.,{x) > 0) 

In example 13. 1[ we have seen that it is possible to produce semi-certificate of 
positivity, where the considered semi-algebraic covering of necessarily have more 
than one piece. 

Remark 3.4. Note that if we restrict ourselves to certificates where the Uij{x) are 
constant matrices, then we a strictly smaller class of certificates. For instance, 

consider the matrix polynomial ( ^ ^9 I which is psd on all R-^, which is obvi- 

\ xy y^ J 

ously psd on a neighbourhood V of the point (1,1) (we may proceed likewise the 
neighbourhood of any given point). 

Let us deshomogenize by setting y — I, set for simplicity V = [l.l + e[, e > 0, 
and assume that 

where the /^'s and gj^s are univariate polynomials psd on V and a; G 
Let 

fi{x) = a, + b,{x-l) + c,{x-lf 
gj{x) = dj + ej{x - 1) + fj{x - if 



Then 



= T.ia1f^(x)+Y.j9J{x) 



Combining these three identities we get 



^(a.-l)V.(:r)+Eff^-(^) = (^-l)' 

i j 

By assumption, we have > and dj > 0. Then, for all j we have dj = and 
Cj ~ 0. Furthermore, if at 7^ 1, then — and also 6j = 0. But substracting the 
last two equalities of the previous system yields 

x-l = ^{ai ~ i)fi{x) 
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a contradiction. 

4. Main results 

There is one case when semi-certificates exist : when the matrix polynomial is 
positive definite. Remind that it means that it is positive definite on the unit sphere 
or the associated projective space. 

Theorem 4.1. Let A{x) G (R[a;])™^™ be a positive definite marix polynomial. Then, 
there is a finite semi- algebraic covering o/R" = U^^iSi, some forms fi.j{x) G R[x], 
some matrix polynomials Ai,j{x) e (R[a;])™^'" such that, 

\fx G S.„ A{x) = 

with the condition that fi,j{x) > for all x Cz Si. 

Proof. Assume that the entries of A{x) are d- forms (d is even). 
We start with a lemma 

Lemma 4.2. Let a{x) be a positive definite d-form and c{x) be a positive semi- 
definite d-form. Then, there exists e > such that a{x) — ec{x) remains positive 
definite. 

Proof. The d-form c{x) = c{x) -\- xf + . . . + xf^ is positive definite and the subset 
— {x \ c{x) = 1} is compact. Since a{x) > on Sc we have for all x G Sc, 
a{x) > m > 0. Hence, for all x, 

Cx X \ 

Then a{x) — ^c{x) is positive definite, which concludes the proof since c{x) > 
c{x). □ 

Let A{x) = {ai,j{xj)i<i,j<m- For a given xq G R" \ {0}, the matrix A{x) is 
positive definite in a semi- algebraic subset Uxg which is open in the unit sphere 
of R". 

Let U be the square matrix of size m whose all entries are equal to 1. Bv l4.2l and 
up to resizing Uxo, we may assume that B{x) — A{x) — €{xf + . . . + xf^)U remains 
positive definite in Uxq for e small enough, and moreover that all entries of B{xq) 
are non-zero. 

Let B{x) = {hj{x)) 

i<i,j<m- Assume that b2.i{xQ) > 0. We may write B{xq) — 
b2,i{xQ)C where C G R™^™ is positive definite. Let hxo{x) be a d-form satisfying 
hx„{xo) = 1. Then, byS^l 

b2.i{x)C - eihxQ{x)Idm 

remains positive definite for some small ei > 0. 
Now write 

B{x) = {b2.i{x)C - eihx„{x)ldm) + B{x) 

Since B{xq) = eildm, the matrix B{x) is positive definite in an open semi- algebraic 
neighbourhood of xq which we still denote by Ux^ ■ 

Then, B{x)) = {bij{x))i<ij<m is such that 62,1(2:) — 0. 
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Of course, we proceed likewise when 62,i(a;o) < 0, writing B{xo) — ~b2,i{xo)C 
where C G R™x™ is positive definite. Thus, 



Then, we repeat the same process to get rid ofi' all the bij{xys such that i ^ 
and reduce, up to resize Uxo, to the case where B{x) is diagonal. Namely, there is 
an open semi-algebraic subset Uxo ^ such that for all x G Uxg, 



where Uk C R™^™ is a constant matrix and fk{x) is a d-form which is psd on Uxg- 

To conclude the proof, we extract a finite semi-algebraic covering by compact- 
edness of S*"^^. 

By homogeneity, we extend the identities ^ to all x G M". Moreover, we may 
manage in order that the polynomials describing the C/^o's are forms, and hence 
the identity ([6]) is true on Sxo, the cone in M" with origin and basis Uxo- Thus, 
we get a finite covering of R" . □ 

Note that the proof gives an open semi-algebraic covering, and that the U^s are 
constant matrices. 

Likewise, we have an analogeous result for biforms : 

Theorem 4.3. Let f{x,y) G M[x,?/] he apositive definite biform of type (n, di; m, ^2). 
Then, there is a finite semi-algebraic covering o/M" = U^^iSi, some forms fij{x) G 
R[a;], some biforms gi,j{x) G R[a:,2/] such that, 



with the condition that fi.j{x) > for all x Si. 

Proof. We proceed as in the proof of Theorem 14.11 We get rid of any monomial 
appearing in / with some odd power with respect to at least one indeterminate. □ 



We say that V is an orthant-neighbourhood of .tq if it contains the intersection of 
a neighbourhood of xq and an orthant Oxo centered at xq (i.e. a subset defined by 
an open condition x ^ xa + (^1, • ■ • , Xn) G Ox„ if ei^i > 0, . . . , e,iX„ > where 
(ei, . . . , e„) G { — 1. +1}"). We will write that Oxo an orthant-neighbourhood of 
(a;o,e) 

The orthant-neighbourhoods fit naturally with the use of Taylor expansion for- 
mula. Let us recall it relatively to an orthant- neigbourhood (xq , e) for a matrix 
polynomial A{x) whose entries are not necessarily homogeneous polynomials : 



B{:x) = {-b2,i{x)C - eihxo{x)Ura) + B{x) 



(6) 




k 



Va; G S^, f{x,y) = ^ ft,j{x){gij{x,y))'^ 



5. Semi-certificates on orthant-neighbourhoods 



^^Ai-\x,) 



a. 
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where we use the standart multi-index symbol for products. Under a condition 
of domination by the constant term, we will see how to derive some certificate on 
orthant-neighbourhoods. 

Definition 5.1. Let A{x) he a matrix polynomial. Denote byV^g the set of all multi- 
indexes P which are minimal (for the lexicographic ordering) such that A'^^^\xq) 7^ 0. 

We say that A{x) satisfies the domination condition at xq if for all multi-index 
a there is some multi-index (3 € T^o o,nd a non negative real number rajj such that 
13 < a and r^^^A^'^) (xo) ± A^"-\xa) is psd. 

Note that T^q — {0} when A{xo) ^ which simplifies the domination condition. 
Note also that if /3 G Txf, then necessarily all its coordinates are even integers. 

Finally, mention that to check positivity via the Taylor formula, it would be 
enough to have the domination condition for any multi-index a whose at least one 
component is odd. 

A typical use of this domination condition appears in the following situation : 

Proposition 5.2. Let A{x) be a psd matrix polynomial on a neighbourhood of xq. 
Then, A(x) satisfies the domination condition at xq if and only if it admits a semi- 
certificate at Xq where the Uij{x) 's appearing in are constant matrices. 

Proof. Assume that A(x) satisfies the domination condition at xq. Let V be an 
orthant-neighbourhood of {xq, e). We have 

where is a subset of all multi-index a such that a > j3 (we have to be carefull 
that one multi- index a may be greater than several elements of r^:^ ) . 

A(xo + eX) = E,.r.„ ((e^^^(^)(xo)) (l - Ea, ^X<^~^r^,,) 

Since = 1, we obtain a semi-certificate on the orthant-neighbourhood V . 
Conversely, assume that 

i 

where each Vi is a constant psd matrix and each fi is non negative on a neighbour- 
hood of Xq. We write the Taylor expansion of fi at xq : 

/,(xo + x)v, = fr/'"^(^")) ^ 

For the matrix polynomial fi{x)Vi, consider the set T^g. If /3 G ^xo: then 
fi^\xo) > and for any a > /3, we clearly have the existence of a positive real 
number r^.^^ such that ra^pf^^\xQ) ± f'f\xa) > 0. This is the domination condi- 
tion. □ 

Remark 5.3. Since positive definite matrices obviously satisfy the domination con- 
dition, we may recover a version of Theorem l4.1l with orthant-neighbourhoods. But, 
maybe (highly heuristic !) it will produce a lot more pieces for the covering, since 
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(again roughly speaking) we need 2" ^ orthant-neighourhoods to recover a usual 
neighbourhood. 

5.1. Semi-certificates relative to a semi-algebraic subset. One may natu- 
rally want to extend the framework of semi-certificates relatively to a basic closed 
semi-algebraic subset S. The problem is that the result given in this section does 
not take into account the equations describing S. We mainly use the underlying 
semi-algebraic set rather than the preordering or the quadratic module generated 
by the equations of S as it is desired for a relative Positivestellensatz. In fact, this 
section concerns more the study of local semi-certificates rather than relative's ones. 

Theorem 5.4. Let A{x) G M[a;]™^'" be a homogeneous symmetric matrix polynomial. 
Assume that A{x) is positive definite for all x Cz S, where S is a closed semi- 
algebraically subset of M" defined by homogeneous polynomials. Then, there is a 
finite semi- algebraic covering of S = Ul^^Si, some forms pij{x) € M.[x], some 
homogeneous matrix polynomials Ai j{x) € R[x]™^'" such that, 

Vx e Si, A(x) = Pi J {x)Af.j {x)Ai^j (x) 
with the condition that pij{x) > for all x ^ Si. 

Proof. We may perform the same proof as in Theorem 14. II □ 

As an example of a case where A{x) is not definite, we consider a general matrix 
polynomial of degree at most 2 in a single (non homogeneous) variable x. 

Example 5.5. Assume that A{x) is psd in the neighbourhood of 0^. 

We assume moreover that ^4(0) is not positive definite otherwise we are done 
by 15.41 or by a domination argument. Up to a base change, we may assume that 

A{Q) = ^ g Q y Then, let us write 

M^)=il nV^Yr ^'V^'il" ^2 A / l + aix + a^a;^ b^x + b2X^\ 

^ ' \ J \ bi Ci J \ b2 C2 J \ biX + b2X^ CiX + C2X^ J 

Since A{x) is psd at 0+, we have ci > 0. 

* If ci > 0, then the domination is satisfied and we get the following certificate 

•^<--' ( ; s ) « ( °'.r ) + --^ ( x ° ) 

where a,l3,ii are positive real numbers chosen such that the constant ma- 
trices of the identity are psd. Note that {1 — fix — ax"^) is obviously positive 
on a neighbourhood of 0+. 

* If ci = 0, then A{x) is psd on 0+ if C2 - 6^ > 0. If C2 - 6? > 0, then we 
have the certificate 

A, ^ N 1 \ / 1 - e bxx \ ( a2 + ae 62 \ 

^(") = ^1-"^")U j + ( M 62 C2-^j 

where e and a are positive real numbers such that the last constant matrix 
in the identify is positive definite. 
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Now, consider the case when ci = and C2 — bl — 0. By positivity of A{x) 
we have aiC2 — 26162 > 0. 

The case when 61 = is trivial since we must have 62 = and the 
certificate fohows. If 61 7^ 0, then we get the condition aibi — 262 > 0. 

First assume that ai6i — 262 > 0. We may write 

0-2 



Aix) 



1 + 



1 




{bix) 



1 + ^x 



ib,x) (1 + ^x) {b,xf 



* In the remaining case when ci = 0, C2 — 6^ = 0, 61 7^ 0, aibi — 262 = 0, 
the positivity of A{x) at 0+ says that 026^ — 62 > 0. The desired certificate 
foUows from the identity given in the previous case. 

This inspection of the most elementary situation leads to conjecture that any 
matrix polynomial in a single variable admits a local certificate of positivity. 

We give a proof of this fact, although we do not give explicit formulas depending 
on the entries as in the previous worked example. 

Theorem 5.6. Let M € R[a:]"^" be a symmetric matrix polynomial whose entries 
are polynomials in a single variable x. Assume that M is psd on a neighbourhood 
ofO+. 

Then, M admits a semi- certificate of positivity at 0^. 

Proof. If A/(0) is invertible, then we are done by 14.11 Hence, from now on, we 
assume that det(M(0)) = 0. 

Let us consider the Smith normal form of AI : 



where E and F are invertible in 



M = EDF 
Jd"^" and D 



diag(cii, . . . ,dn) is diagonal in 



M[a;]"'<" with diagonal entries (di, . . . By changing M to {F^y^MF''^, we 

may assume that Al — ED. Moreover, if = then we may argue by induction on 
the size of M, hence we will assume in the following that > at 0+ and vanishes 
only at 0. 

The matrix M has the form 



M 



( dieia ^261, 2 
^261, 2 ^262, 2 



dn^l.n \ 
dne2,n 



Hence the matrix E has the form 

/ ei,i ei,2 



dn^n^n / 



E = 



'61 2 



62,2 



ei,„ \ 

62.71 



First, let us introduce some notations. 

Define by induction the following sequence of integers 



^ 62, n e^i^ri / 
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Let fco = 1 and set ki to be the first integer j > ki-i such that -r^(0) = 0. It 
yields an increasing sequence of integers : 

1 ~ k^ < ki < . . . < kr < n = kr+i — 1. 

We divide the matrix M = {mk^i)i<k.i<n into block matrices : M = {Mi^j)Q<i,j<r, 
where we set Mij = {rn^kf )i<k<k,+i-k,,i<i<kj+i~kj with rn^^'p = nik+k,-i,i+kj-i- 
Likewise, we divide E — {Eij)o<i,j<r into similar block matrices. 

Let us define the block matrix Mi — {Rk.i)o<k,i<r such that Rk^i — Mk^i when 
(1 < A; < z and I = i) or {1 < I < i and k = i) and Rk,i = otherwise. Namely 

/ ... Ml,, ... \ 

... M,_i,, ; ; 

Mi= ■ ■ . 

M,.i ... Mi,,_i Mi,, : : 

... 

V ... / 

Let us write now 

M = dfc„Mo + dfeiMi + . . . + dk,.Mr. 

Observe then by construction that E{0) is block upper triangular. Since E{0) is 
invertible we get that Ei^i{0) is invertible for each i = 0, . . . , r. 
Moreover, we have 

fc,+i-i 

(7) det(Af,,,) = det{E.u) x [| 

j=ki 

By definition of the /c^'s, it implies that det(^Mi,i) does not vanish at . Hence 

we may derive a kind of domination condition for the matrix M at 0+ . Namely, let 
us write 

M = dk„Mo - dki X fii X Idfeo,fei - . . . - dfe, x x Idk^M 

+dfei(Mi + ^j.ildk„,ki) - dk2 X fi2 X ldki^k2 - . • . - 4-, x A^r x Idfe^^fe^ 



+ dkAMr + firldko,k,) 

= dkoNo + dk,Ni + ... + dk,Nr 

where the /z^'s are positive real numbers and Idp,g is the diagonal matrix in M"^" 
whose i-th diagonal entry is+lifp<i<(7— 1 and otherwise. Because of ([7]), 
we know that we can choose the /i^'s such that the matrix polynomials Ni are psd 
on a neighbourhood of 0+ . 



14 



RONAN QUAREZ 



Then, using the same argument as in Theorem 14.11 we complete the construction 
of a semi-certificate. □ 

The resuh is no more true with more than one variable : 

Proposition 5.7. Let 

1 + ~ y"^ —X 



M{x,y) 



Then, M{x,y) is psd for x^ > y^ > 1 and x^ < y^ < 1 althought it does not admit 
a certificate in a neighbourhood of (!"*", I"*"). 

Proof. First note that det(M(a;, y)) = (?/^ — l)(a;^ — ?/^) to conclude to the positivity 
domain of M{x, y). 

For convenience, let us consider the following change of variables y = 1 + Y 
and X = Y + H . We also perform a base change to "simplify" the expression of 

M(0, 0) = ^ ^^"^ ^ . Let P = ^ I ^ and Q ^ | ^ Then, 



4 




P X Af(0, 0) X Q = 

and thus let us introduce the following matrix 

N{X, Y) = Px M{x, y)xQ 

4: + AH + AY + H'^ + Y^ + 2HY 2H-2Y + H^ -Y^ + 2HY 
2H -2Y + -Y^ + 2HY +Y^ + 2HY 

The matrix N{X,Y) is psd at an orthant-neighbourhood (0+,0+) with respect to 
the variables {H,Y). Let us assume that it has a semi-certificate of positivity. 
Namely : 

where all Fi's and G are psd on an orthant-neighbourhood (0^,0+). The degrees 
of the (2, 2)-entries shows that each fi has only monomials of degree 2 and each Vj 
has monomials of degree 1. Hence : 

F, = a,H^ + b,Y^ + c,HY 
Uj^lj+i^jH + UjY 
Vj = ajH + l3jY 

G = a'H^ + b'Y^ + c'YH + d' H + e'Y + /' 
Identifying the (2. 2)-entries, we get 



(9) 



^c, + 2^a,/?, =2 



Identifying the coefficients in H and Y of the (1, 2)-entries and the constant 
coefficients of the (1, l)-entries yields : 
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(10) 



i 
i 

E7| + /'-4 



Let us introduce a = (oj), /? — (Pj) and 7 = (7^), vectors of the standart 
euclidean space (whose dimension is equal to the number of indexes j). 
After (Uni) and ©, we get 

7^ < 4, < 1 and a • 7 = 2. 

By the Cauchy- Schwartz case of equahty, we get 7^ = 4 and = 1. Thus, 
f ~ and tti = which means that for ah i, at — 0. Moreover, the vectors a 
and 7 must be cohnear and hence a = ^7. 

Likewise /? = — ^7 and — 4. Moreover, for ah i we have 6^ = 0. And hence 

we may assume that there is only one index i and we may write 



(11) N{X,Y)=AHY 



E 



G 



1 




Identifying the (1, 2)-entries, and setting /i — (fij) and 1^ = (vj) we have 



(12) 



( 1 






• M 


1 




2 ' 


• 9 


7 
. 2 ' 





We readily deduce that ^ ~ \- Identifying the (1, 1) -entry, we have furthermore 



' 2 7 • + d' = 4, 
27 • p + e' = 4, 
(13) \ 2fi-P + c' + l = 2, 

p,'^ +a' ^ 1, 
D^+b' = 1, 

By ([T^ and ([T^ . we immediately have d! — e' = 0. Moreover, we get also 

7 • /2 = 2, ^2 < 1 and 7^ = 4. 

By the Cauchy-Schwartz case of equality, we deduce that a' — Q and /I = ^7. 
Likewise 6' = and = ^7. 

Then, G = c'HY where necessarily c' > 0. To conclude, it suffices to note that 
2/2-P + c' + l = 2 gives 2 + c' + 1 = 2, a contradiction. □ 

Remember that any psd biquadratic form in dimension 2 is a sum of squares, 
and hence admits a semi-certificate. This example shows that the local counterpart 
is no more true. 
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In general we may mention the following result. Although completely elementary, 
it looks very much like the ones we can find in [7] or [8] for hermitian squares or 
more classical sums of squares. 

Proposition 5.8. Let A{x) € M[x]™'^'" be a symmetric matrix polynomial. Then, 
A{x) is psd if and only if for all real e > the matrix polynomial 

A, {x) = A{x) + e xf^ Id„ 

admits a semi- certificate. 

Of course the complexity of the semi-certificate may increaes as e goes to zero. 

Remark 5.9. The complexity of a semi-certificate shall be measured by the size of 
the semi-algebraic partitions and also by the number of squares. Here the degrees 
of the polynomials appearing in a (homogeneous) semi-certificate are bounded by 
the degree of the given matrix polynomial. 

Remark 5.10. By [6] with respect to the (compact) unit sphere {X)r=i ~ then 
we obtain for M a Positivestellensatz which is no more homogeneous and where 
the degrees are no more bounded. 

Remark 5.11. Consider a semi-certificate, and a compact semi-algebraic subset Si 
of the given partition. If we assume that the desciption of 5"^ is archimcdean, then 
we may deduce a "true" algebraic certificate without denominators for the matrix 
polynomial on Si, using for instance Schmiidgen or Putinar certificates. 

6. Around the Choi counterexample 

Let us consider the counterexample of a biqudratic psd non sum of squares given 
inn : 

(a;2 -I- —xy —xz \ 

—xy y^ + 2x^ —yz 
—xz —yz z^ + 2y^ J 

Altough it is not a sum of squares, it admits a semi-certificate. Indeed, when 
\x\ > 1^1, it can be decomposed as 

x^ —xy —xz \ / \ 

—xy y^ yz +20 —yz 

—xz yz / V / 

= {-X, y, zY{-x, y, z) + 2(0, -z, ?;)^(0, -z, y) 

/ioo\ /000\ 

-|-2z2 + 2(a;2-z2) o 1 

\000/ \000/ 
And by symmetry, we deduce an analogeous certificate when \z\ > \x\. 

Let us consider a variation around this example ; in the following proposition, 
the Choi counterexample is just Mi : 
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Proposition 6.1. Let 

(x'^ + {X + l)z^ —xy —xz \ 

-xy 2/2_^(A + l)x2 ~yz 

—xz —yz + (A + J 

where A G M. 

Then, M\ is never a sum of squares and it is psd if and only if X (z [0, +oo[. 
Moreover, det(Mo) is a psd non sos polynomial and M\ admits a semi-certificate 
of positivity for any A €]0, +00 [. 

Proof. The exact argument as in ITl works to show that any M\ is not a sum of 
squares. We may reproduce it for convenience. We translate our ramewok into the 
language of biforms, setting 

fj^j^ = {x^+{X+l)z^)S'^-2xyST-2xzSU+{y'^+{X+l)x'^)T^-2yzTU+{z'^+{\+l)y'^)U'^ 

Assume that /ma = ff where fi is a bihnear form. Since there is no monomials 
x^U^, y^S^, z'^T'^ in /maj there is no such monomials in each ff, and hence no 
monomials xU, yS nor zT in each fi. We may write fi = gi + hi were gi depends 
only in the monomials xS, yT, zU and hi depends only in the monomials xT, yU, 
zS. 

Then, /a/^ = Y^{gi + hi)"^ shows that 

^ x^S'^ + y^T'^ + z^U^ - 2{xyST + yzTU + xzSU), 

i 

a contradiction since the right hand side of the equality takes negative values for 
x = y = z = S = T = U = \. 

The fact that Mq is psd and not a sum of squares may also be shown directly by 
using the following consequence of the Cauchy-Binet formula. 

Lemma 6.2. Let M{x) be a symmetric matrix polynomial. If M is a sum of squares, 
then its determinant is also a sum of squares. 0. 

Proof. We first state the Cauchy-Binet formula. 

Given matrices A g E™'*" and B e M-*^™, the Cauchy-Binet formula states that 

det(ylB) = ^ det(^5) det(Bs) 

s 

where S ranges over all the subsets of {l,...s} with m elements, and As (re- 
spectively Bs) denotes the matrix in M™^^'™ whose columns are the columns of A 
(respectively whose rows are the rows of B) with index from S. 

If M(x) is a sum of squares, then it can be written M (x) — A'^(x)A{x) for some 
matrix polynomial A{x) G Then, 

det(M(a;)) = J^s dct{A{x)'^)s det{A{x)s) 
= Es detiMsixf) det(Ms) 
= Es(det(M5))2. 

□ 



The Cauchy-Binet formula shows even more : if M is a sum of squares, then all its principal 
minors are sums of squares. 
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Note that Mq has as a determinant : 

det(Afo) = x^y'^ + y^z^ + z^x^ - ix^y'^z^ 

a variation of the celebrated Motzkin polynomial which is psd and not a sum of 
squares ! 

Thus, the matrix polynomial Mq, which is clearly psd since all its principal 
minors are psd, cannot be a sum of squares since its determinant is not ! 

Thus, Mq (and thus all M\ for A > 0) are psd matrix polynomial. Moreover, 
det(MA(l,l,l)) = \{\ + -if 
which shows that Mx is not psd for small A < and hence for all A < 0. 

Let us study now the existence of semi-certificate of positivity. 

For A > 0, the set of real singular points of det(MA) have projective coordinates 
is [1 : : 0], [0 : 1 : 0] and [0:0:1] (they correspond to the points where the non 
negative polynomial Aet{M\) vanishes). The identity 

( -xy -xz\ [ \z'' \ / \ 

Ma= -xy y2 yz \ + \ Q ^ -2yz \+-{Xx^ ~Az^) \ 10 

\ -xz yz z^ J \ ^2yz Xy^ J \ / 

gives a semi-certificate at the neighbourhood of [1 : : 0]. We may also proceed 
likewise at the neighbourhood of [0 : 1 : 0] and [0:0:1]. 

Elsewhere the matrix Mx is positive definite so that we can argue as in the 
proof of Theorem [4T] and obtain a finite semi-algebraic covering, and hence a semi- 
certificate. 



Whereas, the matrice Mq has a lot more singular points, namely : [1:0:0], 
[0 : 1 : 0], [0 : : 1], [1 : 1 : 1], [1 : 1 : -1], [1 : -1 : 1], [1 : -1 : 1].^ ^ 

So, the problem which remains is to find a certificate of positivity for Mq at the 
points [1:1:1], [1 : 1 : — 1], [1 : — 1 : 1], [1 : — 1 : 1]. It does not appear trivial. 
For instance, the domination condition is not satisfied here. Another natural idea 
would be to consider : 

—xz 

Mq = I —xy y^ + x^ —yz 
















•( 


-xy 







—xz 







) 






y' 






But, setting z = 1 in the last matrix, we get a matrix polynomial which is psd 
at an orthant-neighbourhood of (1+,1+) with respect to the variables {x,y), but 
which does not admit any semi-certificate as shown by Proposition 15.71 

Nevertheless, let us see in the following how it is possible to get a semi-certificate 
in some orthant-neighbourhood of [1 : 1 : 1], and hence a partial result. 
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Let z = 1 and x=l + X,y=l + Y. Write the Taylor formula at (1, 1) : 

Mo = Ao+ XAi + X'^A2 + YBi + ¥"^82 + XYC2 
Consider the following matrices : 



Ax = —Ao + XAi+X^A2 = 



( I + 2X + X2 



"12 



5 Y 



5 Y 



+ 2X + X^ 



_5_ 
12 



__5_ 
12 



AxY = Y^Ao + XYC2 



One may check that 



' 6 



Ay = Y^Ao + YBi + Y^B2 = 



12 ^ 



12 \ 



T2 l+'^y+Y^ 



12 ^ 



\ 12 



12 



1 

3 

-l-XY 
_ 1 

6 



-l-XY 

6 



■n 



= Ax - ^0 



X 



X' 



\ 



6 


12 


12 


-12X^-9Js:-5 


18X^+9X+5 


-12X^+3X-5 


12 


6 


12 


-12X^-9X-5 


-12X^+3X-5 


12X^-3X+5 


12 


12 


6 



is positive semi-definite on all M. Since it is a psd matrix polynomial with respect 
to the single variable X, it is a sum of squares. 
Likewise, one may check that 



V 



12Y^+3r+5 


-12Y^-15Y-5 


-12Y^-3Y-5 


6 


12 


12 


-12Y^-15Y-5 


i8r^+i5y+5 


-12F^-15F-5 


12 


6 


12 


-12Y^-3Y-5 


-12Y^-15Y-5 


18r^+15F+5 


12 


12 


6 



/ 



is positive semi-definite on all M. Since it is a, psd matrix polynomial with respect 
to the single variable F, it is a sum of squares. 
Finally, write 

'2{1 + 6XY) 



AxY = 



12 



Ao - XY{Ao - C2) 



Where wc check that (^o — C2) is a psd constant matrix. 

If we sum up all these informations, wc obtain a desired certificate of positivity 



Mo=Ax+Ay + Axy = Aq 



X 



-X^]+Cx+A 



-Y 



-Y^]+Cy + A 



XY 
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with respect to the open orthant-neighbourhood of (0,0) defined hy X > and 
y < 0. 

Likewise, we may produce a local certificate with respect to the orthant defined 
by X < and y > 0). 

But, it seems less clear how to obtain a local certificate with respect to the 
orthants defined hy XY > 0. 

□ 

7. Concluding remarks 

7.1. Open questions about semi-certificates. We have introduced the notion 
of piecewise semi-certificate of positivity for matrix polynomials. For the moment, 
the only general result is that all positive definite matrix polynomial admit such a 
certificate. A lot of things remain to be studied. 

• We shall better understand the set of all psd matrix polynomials which 
admit a semi-certifiacte. Beginning with biquadratic forms for instance ? 

• We shall developp some effective algorithm to produce the certificates. 

• About the complexity of certificates : how can we bound the number of 
squares, and also the number of pieces of the semi-algebraic partition ? 

On the other hand, one may also rise some questions about the familly of psd poly- 
nomials which have a quadratic determinantal representation. This is the object of 
the last subsection : 

7.2. Semi-definite quadratic determinatal representations. In the spirit of 
what happens in Proposition 16. 1[ one may be interested in determining what psd 
polynomial can be written as the determinant of a psd quadratic matrix polynomial 
(whose entries are homogeneous polynomials of degree 2). Namely, 

Question 7.1. Let f{x) G M[x] be a form of degree 2d which is suppose to be non 
negative. Does-there exists a quadratic matrix polynomial M{x) G R[x]'^'^''^ which 
is psd for all x G R" and such that 

(14) f{x) = det(Af(a;))? 

One motivation for such a question is that it appears as the quadratic analo- 
geaous of very classical linear determinantal representations which have been stud- 
ied for a very long time. The specifically real considerations being more recents and 
esentially du to Vinnikov (see for instance [H] and all related papers). Let us re- 
call that these have a lot of applications for instance to Linear Matrix Inequalities, 
convex modelling, etc... 

On the other hand, biquadratic forms are a quite common object found at various 
areas in ingenering applications. For instance, any determinantal representation as 
in (|14p of a psd non sum of squares polynomial provides, by Cauchy-Binet formula, 
another example of a psd biquadratic form which is not a sum of squares. 

Of course, a simple count on the number of parameters shows Question 17.11 has 
a negative answer in general. 

We can be even more precise : if / is a form of degree 4 in 4 variables, such 
a representation never exists when / is psd non sum of squares. Indeed, if / = 
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det(M(a;)) with AI{x) G ]R[a;]^^^, then M{x) is a psd biquadratic form which is a 
sum of squares (in dimension 2), hence det(M) is a sum of squares of polynomials 
by Cauchy-Binet formula. 

We may also look at the trivial case when n = 2, which we may deshomogenize 
for simplicity. The answer to Question 17.11 is clearly positive by an elementary 
argument. Indeed, let us write 

2d 

fix) = = - "^■)' ^ + ^'^)' + ^fe) 

i—O j k 

Then, the result follows obviously from the multiplicative property of the determi- 
nant and the fact that the (x — ai)^'s and the (a; + + 7^)'s are psd quadratic 
polynomials. 

We may even give a more algorithmic construction (a polynomial time algorithm 
with respect to the size of the coefhcients of the polynomial) using arrows matrices 
and following [5] and [9]. 

Back to Question 17. 1[ and before expecting general results, one may wonder if 
for instance the celebrated Motzkin and Robinson polynomials 

{Mo = + x'^y'^ + x^y"^ - 3 x'^y'^z^ 
RO = + y6 + ^6 _ ^4y2 „ ^4^2 _ y4^4 „ ^4^2 „ ^4^2 „ ^4^2 ^ 3 ^2y2 ^2 

can be written as in ([T4|) ? 

It is quite easy, using for instance the linear well-know determinantal representa- 
tions of cubics curves, to produce quadratic determinantal representations for these 
two polynomials, but unfortunately none such representation yields a psd quadratic 
matrix polynomial. 
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